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Abstract 

The discrete-time version of totally asymmetric simple-exclusion process (TASEP) on a finite 
one-dimensional lattice is studied with the periodic boundary condition. Each particle at a site 
hops to the next site with probability < p < 1, if the next site is empty. This condition can be 
rephrased by the condition that the number n of vacant sites between the particle and the next 
particle is positive. Then the average velocity is given by a product of the hopping probability p 
and the probability that n > 1. By mapping the TASEP to another driven diffusive system called 
the zero-range process, it is proved that the distribution function of vacant sites in the stationary 
state is exactly given by a factorized form. We define /c-particle velocity correlation function as the 
expectation value of a product of velocities of k particles in the stationary distribution. It is shown 
that it does not depend on positions of k particles on a circle but depends only on the number k. 
We give explicit expressions for all velocity correlation functions using the Gauss hypergeometric 
functions. Covariance of velocities of two particles is studied in detail and we show that velocities 
become independent asymptotically in the thermodynamic limit. 
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I. INTRODUCTION 



The one-dimensional totally asymmetric simple- exclusion process (TASEP) is a minimal 
statistical-mechanics model for driven diffusive systems of many particles with hardcore 
exclusive interaction jl-6|. In the present paper, we consider a discrete-time version of the 
TASEP on a circle, i.e. a one-dimensional finite lattice with periodic boundary condition, 
in which the parallel update rule is applied This process can be exactly mapped to 
another driven diffusive particle system called the zero-range process (ZRP), by regarding 
each number of vacant sites between successive particles in the TASEP as a number of 
particles at each site in the ZRP js, 9]. 

In the ZRP representation, particles hop from site to site on a lattice with a hopping 
probability which depends only on the number of particles at the departure site. We assume 
that there do not occur any creation and annihilation of particles in the ZRP. Then, if a site 
is vacant, there is no possibility that a particle hops from that empty site to other sites, as 
a matter of course. The configuration in the TASEP such that the next site of a particle is 
occupied by another particle is represented by the configuration in the ZRP such that the 
corresponding site is empty. Therefore the prohibition of hopping by hardcore exclusion in 
such a jamming situation in the TASEP is automatically satisfied in the ZRP. Moreover, 



the steady state^o 
factorized form 



' the ZRP is exactly described by the probability density function in a 
9| , and then the stationary distribution of vacant sites in the TASEP on 
a circle is explicitly determined as Eq.flH]) given below js il]. 

The map from the TASEP to the ZRP is interesting from the viewpoint of quantum 
statistical mechanics, since it seems to be a map from an interacting Fermi gas to a free 
Bose gas. On the other hand, as demonstrated in the present paper, the procedure in which 
we analyze the TASEP/ZRP is quite different from the standard way for free Boson systems 
in the following sense, (i) Thermal equilibrium state of free Bose gas is usually treated in 
the grand canonical ensemble by introducing fugacity, while here we want to consider the 
driven diffusive system with a fixed number of particles and thus we treat the system in the 
canonical ensemble. Then the canonical partition function Z^, where L and iV denote the 
numbers of particles and vacancies in the TASEP, respectively, plays an important role in 
calculation, (ii) In the usual theory of free Bose gases, condensation of particles in a specified 
energy level is carefully studied {e.g., the Bose-Einstein condensation in a ground state). In 
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the context of ZRP, however, distribution of vacancies (i.e., empty energy levels) should be 
well studied. The reason is that, if a site in the ZRP is occupied by one or more than one 
particle, then the velocity of the corresponding particle in the TASEP can be positive, but 
if the site is empty in the ZRP, the velocity of the particle in the TASEP is definitely zero. 
When we consider the TASEP as a simple model of traffic flow, the flux which is defined as 
the product of velocity and particle-density is the most important quantity. The purpose of 
the present paper is to study velocity correlations of particles in the stationary state of the 
TASEP on a circle. 

Let N = {1, 2, 3, . . . } and N = {0} U N = {0, 1,2,...}. For K G N, we consider a one- 
dimensional lattice A = {1, 2, ... , K}. Each site i G A is either occupied by a particle, which 
is denoted by r](i) = 1, or vacant denoted by i](i) = 0. The following discrete-time stochastic 
process is considered for simulating the TASEP on a circle. Let < p < 1. At each time 
t G N , given a particle configuration rj t = {r] t (i)}ieA G {0, 1} A , let A t = {(i,i + 1) : 1 < 
% < K, s.t. rj t (i) = 1, rj t (i + 1) = 0}, where the periodicity r](i + K) = r](i),i G A is assumed 
and the nearest-neighbor pair of sites (K,K + 1) is identified with (K, 1). Every particle 
at site % such that [i,i + 1) G A t has chance to move to its next site i + 1, since the site 
i + 1 is vacant; i] t (i + 1) = 0. But, in general, only a part of such particles move depending 
on the parameter p as follows. We choose a subset of A t randomly, in the sense that each 
pair of nearest-neighbor sites + G A t is chosen independently with probability p. The 
obtained subset of A t is written as A t . In the present paper, the total number of elements 
included in a set B is denoted by \B\ and, for B C C, the complementary set of B in the 
set C is expressed by C\B. (By definition \C\B\ = \C\ — \B\.) Then the probability that 
A t is chosen from A t is given by p' At '(l — p)l A *H A *l. Only the particles at sites {i} such that 
(i,i + 1) G A t indeed move to their next sites. That is, the particle configuration at time 
t + 1, rj t+1 = {rj t+1 (i)} ieA , is given by 



Here note that by definition of A t , if {i, i + 1) G A t , then (i — A t . The parameter p is 

called the hopping probability and the above procedure is said to be the parallel update rule. 
The total number of particles is conserved in the process, which we write L in the present 




(1) 
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paper. We assume 1 < L < K. 

Given a particle configuration r\ G {0, 1} A , let %\ = min{l < i < K : 77(f) = 1} and define 



= min{ij < i < K : = 1}, 1 < j < L — 1, 
i.e., ij is the site occupied by the j-th particle, 1 < j < L. Then we put 



n{j) = ij+i - ij 



1<J<L, 



(2) 



(3) 



where ii+i = h+K. That is, n(j) gives the number of vacant sites between the j-th and the 
(j + l)-th particles. By Q with ([2]), a configuration of vacancies n = {n(j)}j =1 is uniquely 
determined from the particle configuration rj = {f]{i)}f =1 . 

We should note that n does not determine rj uniquely, however, since the information on 
the position of the first particle, i±, is missing in the map 77 — > n. This information may 
be, however, not important, since here we consider the TASEP on a circle. The stochastic 
process nt=_\n t (j)}f =1 , t G N , obtained from r] t ,t e N by this map, is a special case of 
the ZRP 8|, l9( . As a consequence of general theory of ZRP the probability distri- 

bution function in the stationary state Pl i at(ii) of configuration n of vacancies is uniquely 
determined as follows. Since the lattice size K and the total number of particles L are 
conserved, the total number of vacant sites iV = ^f=i = K — L is also a constant. We 
fix 1 < L, N < K. Then the configuration space of n is given by 

n L , N ={n = {n(j)}f =1 G {0, 1, . . . , N} L : ]T n(j) = N\ (4) 



and 



where 



L,N{n) 



13| 



1 L 

_J]/(nO-)), n = {n(j)}f =1 en L!N , 

1L > N 7=i 



/(") 

and the partition function is given by 

Zl,n = 



1^) 



n-l 



if n = 0, 
if n > 1, 



(5) 



(6) 
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e n 



l-p^LF ( 1 - JV;2; 



(1 -p) i+r 



F ( L + 1, AT + 1; 2; 



I — p 
1 



(7) 



with the Gauss hypergeometric function 



F( a ,/3; 7 ;.) = £M^^, (8) 
a) = 1, (a) n = a(a + !)••• a + n — 1), n > 1. Note that the last equality in ([7]) is due to 



Rummer's transformation 



F(a, (3; T , *) = (1 - zf-^F^ - a , 7 - # 75 

We write the expectation with respect to the stationary distribution ([5]) as E^^f-] in this 
paper. 

For the j-th particle, 1 < j < L, if n(j) > 1, that is, if n(j) ^ 0, then that particle can 
move to the next site with probability p in a time-step. Then if the velocity of j-th particle 
is denoted by Vj, the expectation of this random variable in the stationary distribution P^at 
is given by 

E LiN [Vj] = pE L , N [l(n(j) > 1)] 

= P E 1 ( n 0') ^ X ) p ^( n )' ( 9 ) 

where l(w) is an indicator of an event u; l(u>) = 1 if u) occurs, = otherwise. The 
average velocity (jUJ) is independent of j, since the system is homogeneous in space, and it 



has been explicitly calculated as |10|, 111 ] 



EN-l/i \N+l-n y 
, , n=0 \ L ) Z/ ^." 

-^L.TVlVJ - 

P F(l-L,l-iV;l;l/(l-p)) 
LF(l-L,l-N;2;l/(l-p)) 
1-p F(L,N;l;l/(l-p)) 
L F(L+l,N+l;2;l/(l-p)Y 

The density of particles is given by 



(10) 



and the flux Jl^ is defined by 

J L , N = P E L>N [V}. (12) 



If we plot Jl,n versus p, we obtain a fundamental diagram as demonstrated by Kanai [11] 
(see Fig{T]in the present paper). Moreover, Kanai et al. [k]] determined the thermodynamic 



Flux 
0.06 



0.04 



0.02 




FIG. 1: The fundamental diagrams. The flux Jl,n is plotted versus p for K = 20 with p = 0.2 
by a solid line. The broken line shows the thermodynamic limit (|13p with p = 0.2. We find that 
the difference between the results for finite systems and the thermodynamic limit is very small, in 
particular, when p ~ 0, p ~ 1.0. 

limit, i.e. the scaling limit of L — > oo, iV — > oo with keeping p = L/(L + N) be a constant, 
for the average velocity and obtained (see FigJT]) 



lim E^. WhMli!) , <*„<!. (13) 

L-^oo,N^od; Zp 
p=const. 

This result coincides with the exact solution for an infinite system obtained by Schadschnei- 

n 

der and Schreckenberg [151 ] . 

In the present paper, we study velocity correlation functions and show that, as exten- 
sions of the formula (fTDI) . they are generally expressed by using the Gauss hypergeometric 
functions. In particular, the covariance of velocities V and V of two particles at different 
sites, 

Cov L , N [V, V'\ = E LjN [W] - E iiiV [Vf , (14) 
is studied in detail and it is shown that 



lim Gov LN [V,V'} = 0. 

L— >oo,N— >oo; 
p=const. 



(15) 



That is, velocities of particles are correlated in the stationary state in any finite systems, 
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but it is proved that they become independent asymptotically in the thermodynamic limit 
in the discrete-time TASEP on a circle. 

The paper is organized as follows. In Sec. II. A velocity correlation functions are defined 
and the general formula is derived. As special cases, the obtained expressions of average 
velocity and covariance of velocities of two particles are studied in detail in Sec. II. B and 
C, respectively. Sec. Ill is devoted to proving asymptotic independence of velocities ( !T5|) in 
the thermodynamic limit. Concluding remarks are given in Sec. IV. Appendix A is given for 
showing the derivation of a Riccati equation, which governs the covariance of velocities and 
is used in the proof in Sec. III. 



II. VELOCITY CORRELATION FUNCTIONS 



A. General Formula 



Let 1 < k < L. From the L particles, we pick up k distinct particles arbitrarily; let 
1 < jx < J2 < ■ ■ ■ < jk <■ L, and consider a set of k particles such that the £-th particle in 
this set, 1 < £ < k, is originally the je-th particle in the whole particle systems. Write the 
velocity of the ji-ih particle by Vj v 1 < £ < k. The velocity correlation function for the k 
particle is then defined by 

~ k 

El, W l[V Jt 



FL,N{{V je } i=1 



,1=1 



P 



E \{^n{ 3i )>l)V L ^n) 



(16) 



nefl L N i=\ 



Since the stationary distribution function is given by the factorized form (Q, it is written as 



N N 



N k 



FL,N{{Vjg} l=1 

/JL,N n(l)=0n(2)=0 n(L)=0 l=\ 
, N N 



E E •■■ E H 1 ^)^ 1 )!!^)) 1 E^') = iV 



J L,N 



n e /«?)) n e /("(o)i e 



mi) 



N 



jel L \J k n(j)=0 



.1=1 



p 



J L,N 



i&Jk n(£)=l 

n e fw)) n ( e - /(°) 1 1 ( e 



N 



n[i\ 



N 



(17) 



jel L \J fe n(i)=o teJ fe \n(e)=0 

where I L = {1,2, . . . , L} and J k = j 2 , . . . ,j k }- 



8=1 



We perform the binomial expansion 

(N \ N 

E /(»(*)) - /(°) = E (-/(°)) |JAK| n E ( 18 ) 
n(£)=0 J KcJfe <£Kn(()=0 

where the first summation in the RHS is taken over all subsets K of Jfc and | \ K| is the 
number of elements in the complementary set of K in Note that I £ = (I £ \Jfc) UKU (Jfc \ 
K) and i G \ K implies n(i) = 0, since the weight /(O)/^ ^ is assigned in the expansion. 
Therefore we can set J2i=i n (i) = J2ie(i L \J k )uK n (^) * n Since we set f(0) = 1 as ()6]), 

we obtain 

k N N 

P 



E(- i ) |jaki n E^^^riE^^)) 1 e n ^= N 



L,N KeJ fc jei L \J k n(i)=o teKn(«)=o \ie(i£\J fc )uK 



A' 



iV 



^(-i)Wi jj ^; /(„(<))! e ^) = Ar • ( 19 ) 

L,N KGJ fc ie(I L \J fc )UK n(i)=0 \ie(Ii\J fc )UK / 

When |J fe \ K| = s, < s < fc, |K| = k - s and thus |(I £ \ J fc ) U K| = |I L \ J fc | + |K| = 
(L — k) + (k — s) = L — s. Therefore 

n e ( e 

ie(I L \J t )UKn(i)=0 \ie(I L \Jfc)UK) 
L-s N /L-s 

=n e /("(p)) 1 E^) =Ar 

p=l n(p)=0 \p=l 
L-s 

= e n^ n (p)) =z w 

Si.ce the numbe, of distact subsets K in J 4 satisfying | J, \ K| = . is Q , < . < *, 09 

is equal to {p k /Z L N } XL=o( — ^-T Cs) ^l~ s ,n- The result does not depend on the choice of 
particle positions Jfc = {ji, j 2 , ■ ■ ■ , Jfc}, but depends only on the total number k of particles, 
whose velocity correlation is calculated. This special property comes from the factorized 
form (jSJ) of the stationary distribution in the present system, in which the factor f(n) are 
independent of the system sizes, L and N, as given by ([6]). We summarize the result by the 



S 



following formula, 



F L , N (k) = F L , N ({V je } k i=1 



P 



s=0 



J L-s,N 



p 



LF(l-L,l-JV;2;l/(l-p)) 



x 



J2(-1Y(L -s)(^f(i-L + sA-N; 2; ^ 



LF(L + l,iV + 1; 2; 1/(1 — p)) 

k 



(20) 



S J \ v J \ 1 — p , 

It should be noted that still velocities are correlated in the sense that Fi^(k) ^ 
(F LtN (l)) k , k > 2. In other words, 



E 



L,N 



i=i 



^l[E L , N [V h ), k>2. 



B. Average Velocity 



By setting k = 1 in the general formula (12Uj) . we obtain 



E L>JV [y] = F X)iV (i) 



J L,N 



Z L -i. 



A' J 



LF(1 - L, 1 - N; 2; 1/(1 - p)) - (L - 1)F(2 - L, 1 - N; 2; 1/(1 - p)) 
LF(l-L,l-N;2;l/(l-p)) 
LpF(L + 1, JV + 1; 2; 1/(1 - p)) + (L — 1)(1 — p)F(L, iV + 1; 2; 1/(1 - p)) 



LF(L + 1,N + 1;2; 1/(1 — p)) ' (21) 

Now we show that the last expression of (12TT) is equal to the last expression of ( ITOj) . First 
we rewrite the numerator of fl2TT) as 



-(1-p) 



LF L + 1,JV+1;2; 



1 



-(L-l)F L,iV + l;2 



1 — pj 1 — p 
1 



L F f L + 1, AT + 1; 2 1 



1 — p 



1 — p 



(22) 



If we use the recurrence relation of the Gauss hypergeometric series 14 1 

azF(a + 1, p + 1; 7 + 1; z) = 7(^(0, + 1; 7; z) - F(a, (3; 7; z) 
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for the second term in ([22]), (1221) becomes 



-(1-p) 



LF L+l.N + 1:2: 



1 



1 — p 



F [ L,N + !',!', 



1 



+F L,N:1: 



1 — p 



(L-1)F[L,N + 1;2; 



1 — p 
1 



1 — p 



Next we apply the formula 

( 7 - a - l)F(a, /3; 7; z) + uF(a + 1, /3; 7; z) = (7 - l)F(a, #7-1; *) 
to the first and the fourth terms in (1231) . Then the sum of them becomes 



LF L + l,jV + l;2; 



(I-l)F L,JV + 1;2; 



1 



(23) 



F L,JV+1;1; 



1 



-| / \ / I 7 '77-1 I I 7 ' '771 

1-P/ \ 1— £>/ \ \—P t 

which is cancelled by the second term in (1231) . Therefore, the numerator of (12 lj) is equal to 
— (1 —p)F(L, N; 1; 1/(1 — p)), and the equivalence of the last expression of ( 12TI) and the last 
expression of (TX0T) is confirmed. 



C. Covariance of Velocity 



By the first expression in (1201 for k = 2, we obtain 



Vl,n[W] = F LyN (2) 

~2 



V 



[Zl,n — 2Zl-i.n + 



N 



J L,N 



Then the covariance of velocities ( 1T4|) is given by 

Cov LiJV [v;y'] =f 



2 I Zl-2,n ( Zl-i^n 



(24) 



(25) 



Zl,N V ^L,iV 

By the expression of partition function ([7]) using the hypergeometric function, it is written 



as 



Cov LjN [V,V] 



P 



(L-2)F(3-L,l-iV;2;l/(l-p)) 
LF(l-L,l-iV;2;l/(l-p)) 



(L-l)F(2-L,l-iV;2;l/(l-p)) 
LF(l-L,l-JV;2;l/(l-p)) 

(L-2)F(L-l,iV+l;2;l/(l-p)) 
LF(L + 1, iV + 1; 2; 1/(1 — p)) 

(L-l)F(L,iV + l;2;l/(l-p)) 
LF(L + l,iV+l;2;l/(l-p)) 



(26) 
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Cov 




FIG. 2: The exact solution of the covariance of velocities (|26|) is plotted versus the particle density 
p for K = 100 with p = 0.5. Computer simulations are also performed and numerical results are 
dotted. 



If we use the recurrence relation of the Gauss hypergeometric function 



a(l - z)F(a + 1, /3; 7; z) + 7 - 2a + (a - (3)z F(a, (3; 7; z) 
+ (a-7)F(a-l,/3;7;z) = 0, 



the second expression of ( |26|) is rewritten as 

Cov l mv,V] = (1-p) 2 1 



1 — p 

YlAp) 



1+1- 



L-N-l 1 
2(L-1) 1-p 

L-N - 1 1 



1 - 



where 



2(L-1) l-p 
(L-l)F(L,N+l;2;l/(l-p)) 



(27) 



(28) 



LF(L + l,iV + l;2;l/(l-p)) " 
Figure [2] shows Covl,n[V,V] given by ( 1261) for K = 100 with p = 0.5 as a function of 
the density of particles p. This figure shows that velocities are negatively correlated. We 
performed computer simulations for a variety of systems with K = 100 and p = 0.5 by 
changing the number of particles 1 < L < K. For each trial, we discarded first 2 x 10 3 
time-step data and used 8 x 10 3 time-step data for evaluating the covariance of velocities. 
In Figj2j each dot indicates the averaged value of 10 4 trials. We can find that the simulation 
data coincide very well with the exact solution (|26l) . 
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III. ASYMPTOTIC INDEPENDENCE OF VELOCITIES 



A. Riccati Equations 



Here we consider the average velocity ( 1101) as a function of the hopping probability p and 
express it by 

v LtN (p)=E LiN [V}. (29) 
Kanai et al. [lO] proved that it solves a Riccati equation 



dv L , N (p) 



+ 



. , 2 L + N . . N 
-vl,n{p) jz \V L ,n{P) + 



0. 



dp ' p(l — p) ' p{l — p) yi ' ' 1 — p 

From this equation, the thermodynamic limit ( JT3]) was concluded. (See also (|36|) below.) 
We have found that Y L;N (p) defined by ff28l) also solves a Riccati equation in the form 
dY LjN (p) L — 1 Vl-i,n(p) 



(30) 



dp 
+ 



+ 



P v L>N (p) 



-YlApY 



vl-i,n(p) ] 2{L-l)p-{L + N -I) 
dp v LyN (p) p(l-p) 



YlAp) 



L-l v L , N (p) 



0. 



1 -pv L ^ hN (p) 

The derivation is given in Appendix A. Note that by (|30|) we obtain the equation 

d_ j vl-i,n(p) _ dv L -i )N (p)/dp _ dv LtN (p)/dp 
dp vl,n(p) v L - 1)N (p) v l,n(p) 

(L - l)v L -i, N {p) - Lv L , N (p) + 1 + Np 



(31) 



p(l-p) 



• (32) 



Vl-xAp) vl,n(p) J _ 

That is, the Riccati equation f l3~T|) for Y L ^ N (p), which governs Covl,n[V, V] through fT27|) . is 
coupled with the Riccati equations (1301) for vl^{p) and vl-i,n(p)- 



B. Large-size expansion and thermodynamic limit 



Following the procedure given by 10J, we consider the power expansion of the quantities 
with respect to the inverse of system size, 1/K with K = L + N, 



11 

Vl,n{P) =v + v x — + v 2 ^ H , 

v L -i,n(p) = v' + + v' 2 -^ + ■ ■ 
Y L>N (p)=Y + Y 1 ± + Y 2 ±- 2 + ... 



(33) 
(34) 
(35) 
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where the coefficients Vi,v' i7 Yi, 2 = 0,1,2,... are assumed to be functions of p and p. Putting 
fl33|) and (jMj) into (130]) and its modification obtained by setting L — > L — 1, and taking the 
thermodynamic limit if — » oo with p = L/K = const., the first terms in (I3"3"|) and ( I34p are 
determined as [lo| 



w = v 



1 - y/1 - 4 PP (1 ~ P ) 

2p 



(36) 



This result implies (|T3|) . 

Similarly, we put ( I3"3"|) -f l3l)|) into f[3"Tj) with ( 13"2"j) . In the thermodynamic limit, the differ- 
ential equation is reduced to the algebraic equation 



p(l-p)Y 2 -(l-2pp)Y -pp = 



for Y , which is solved as 



1 - 2 PP ) ± y/l - Appjl - p) 

2p(l-p) 

On the other hand, in the similar way, we can show that (j27p gives 



Yq = lim Y LtN (p) 

K — >co; 
p=const. 



(37) 



lim Cov LiJV [V;y'] =p(l-p) + 

if— >-oo; 
p=const. 



1 - 2pj9 



2p 



;i-p)^o 



1 - 2pp 
2p 



(38) 



If we apply the result (I3TI) . turns to be 



lim Cov^jv^V'] = p(l-p) + 

p= const. 

= 0. 



1 - 2pp 
2p . 



1 - 4pp(l - p) 
Ap 2 



(39) 



Vanishing of the covariance implies that velocities of particles of the discrete-time TASEP 
become independent asymptotically in the thermodynamic limit in the stationary state on 
a circle. 

Figure [3] shows the numerical data demonstrating Cov^jvfV, V] — > as K = L + N — > oo. 
Here we set p = 0.5 and p = 0.6 and performed computer simulations by increasing K from 
100 to 1000 by 150, in which Cov L>A r[V, V] with L = pK and N = (1 - p)K are evaluated. 
The linear fitting of Covl,n\Yi V] versus 1/K gives 

rn _ C (P,P) 



K 



(40) 



with c(p = 0.5, p = 0.6) ~-4x 10" 
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FIG. 3: Numerical fitting of Covl^[V, V'] evaluated by computer simulations versus l/K for 
p = 0.5, p = 0.6. The data shows that Covl,n[V, V] becomes zero as K = L + N — >• oo in the form 
(|40p with c(p = 0.5, p = 0.6) ~ —4 x 10~ 2 . In each trial of computer simulation, we discarded first 
2 x 10 4 time-step data and used 10 4 time-step data for evaluating the covariance of velocities and 
each dot in this figure indicates the averaged value over 10 4 trials. 

IV. CONCLUDING REMARKS 

In the present paper, we study a discrete-time version of the TASEP on a circle developed 
by the parallel update rule. This version is easily simulated by a computer. As a matter of 
fact, we have checked the validity of our exact solutions for finite-size systems by comparing 
them with the numerical simulation data as shown by FigO 

From the viewpoint of statistical mechanics, it is important to discuss the thermodynamic 
limit for non-equilibrium steady states realized in the present model. It is rather difficult, 
however, to define parallel update dynamics for a system with an infinite number of particles, 
since the number of updated sites can be infinity, i.e., \A t \ = oo using the notation in 
Sec. I. Kanai et al. [k]] overcame this difficulty and determined the thermodynamic limit of 
average velocity (fl~3l) . They obtained the differential equation which governs the average 
velocity and took the thermodynamic limit in the equation. In the present paper, we extend 
their procedure for the covariance of two-particle velocities and showed that it becomes 
asymptotically zero in the thermodynamic limit (}39l) . 

It should be remarked that both equations which govern the average velocity and the 
covariance are given by the Riccati-type differential equations with respect to the hopping 
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probability p. As a matter of course, the Riccati equation (T3TT) for Yl,n(p), which governs 
Cov^jvfV", V] through ([2~Tjh is coupled with the equations f[30l for the average velocities (see 
(13"2"|) ). and thus it becomes much complicated. Further study of hierarchy in the coupled 
system of differential equations which determine the third and higher-order moments of 
velocities will be an interesting future problem. 
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Appendix A: Derivation of Riccati equation (1311) 



By definitions ([28]) and ([29} with (fTOj). 

Vl-i,n(p) 



vl,n(p) 



Yl,n(p) 



F(L-l,N;l,z) 
F(L,N;l,z) ' 



where we put 



1 



1 — p 



By the following formula of the Gauss hypergeometric function 

'1 _ Z )L+N-1 



d 

dz 



yL-l 



-F(L,N;l;z 



(I _ z) L+N - 2 
1 - L)- L F (L - 1, N; 1; z), 



OAip is written as 



with 



Vl-i,n{p)xt /\ z(l-z) 1 dw(z) 

t^-il,n\P) - — j ~r~\ j 

Vl,n{P) I — L w{z) dz 

{\ - z) L+N ~ l 



w(z) 



yL-l 



-F(L,N; 1; z). 



Here we consider the generalized hypergeometric differential equation 



du{z) y—^ 1 — Aj — du\z) 



dz 2 



E 



i=i 



z — a, dz 



n 

i=l l<j<3;i^i 



u(z) 



(z - a x ){z - a 2 )(z - a 3 ) 



(Al) 



(A2) 



(A3) 



(A4) 



(A5) 
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where Fuchs' relation 5^f =1 (Aj + AQ = 1 is assumed to be satisfied. The solution of (1A5j) is 
expressed by 

(X\ (X2 ^3 

Ai A2 A3 z 
X[ A 2 A3 



u(z) = P< 



(A6) 



which is called Riemann's P function 14j. As a special case, the Gauss hypergeometric 
function (jSJ) is given by 



F(a,/3; 7 ^) = P< 



In general, the following relations holds, 



1 00 

a z 

1 — 77 — a — (3 (3 



P < 



(z- ai ) k P < 



By (|X7jl . (jMD is written as 



Qi G 2 O3 

Ai A2 A3 2 

A'i A 2 A 3 

ai a 2 00 

Ai A2 A3 z 

\[ A 2 A 3 



> = P< 



> = P< 



0,2 CL\ CI3 

A2 Ai A3 z 
A 2 Ai A 3 



ai a 2 00 
Ai + k A 2 A 3 — k z 
Ai + fc A 2 A 3 - 



{z) = {-l) L+N -\z-l) 



\L+N-1JI-L 



1 00 
L z 
1-L-JV iV 



Using the relations (1A8I) and (1A9I) . we can show that 



1 00 

1 — L L + N-l L-N z 
l-L 



It implies that w(z) solves the differential equation 

d 2 w(z) (L-N + l)z-2L + ldw(z) (L — l) 2 



dz 2 



(l-z)z 



+ 



dz (1 — z)z 2 



w{z) = 0. 



(A7) 



(A8) 



(A9) 



(A10) 



(All) 
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Since dw(z)/dz is related with Y^^ip) by (1A3j) . (lAllf) gives the first-order differential equa- 
tion for Yl,jv(p)- By straightforward calculation, fl31~]) is derived. 
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